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Abstract 

A Bianchi type-I cosmological model in the presence of a magnetic 
flux along a cosmological string is considered. The first objective 
of this study is to investigate Einstein equations using a tractable 
assumption usually accepted in the literature. Quantum effects of the 
present cosmological model are examined in the framework of loop 
quantum cosmology. Finally we draw a parallel between the classical 
and quantum approaches. 
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1 Introduction 



The Bianchi models, which describe homogeneous but anisotropic space- 
times, have been extensively discussed in the literature, motivated in part 
by attempts to explain small but significant anisotropics in the cosmic mi- 
crowave background (CMB) [HE] and large structures [3]. 

Among the anisotropic Bianchi models, the simplest ones are Bianchi 
type-I (BI) models whose spatial sections are fiat but the expansion or con- 
traction rates are direction dependent. 

Primordial magnetic fields can have a significant impact on the CMB 
anisotropy. Also the presence of strong magnetic fields raises interesting 
problems like the formation of galaxies in the Universe. The BI models are 
appropriate for the investigation of a Universe which is permeated by a large 
scale, homogeneous magnetic field. 

In the early stages of the evolution of the Universe it is expected that topo- 
logical defects could have formed naturally during phase transitions followed 
by spontaneous broken symmetries. Cosmic strings are linear topological 
defects, have very interesting properties and might play an important role in 
structure formation [U [5] . 

In the first part of the paper we shall investigate the evolution of a BI 
model in presence of a cloud of strings and magnetic field. In order to solve 
Einstein equations we resort to a tractable assumption concerning a relation 
between the rest energy and tension density of the system of strings [6]. 

Loop Quantum Gravity (LQG) [H [8] represents one of the most com- 
pelling attempt towards a complete non perturbative quantum theory for 
the gravitational interaction. The cosmological application of LQG was de- 
veloped in terms of invariant connections [D] and this model was denoted 
by Loop Quantum Cosmology (LQC). LQC takes the ingredients of LQG 
and applies them to expanding Universe or black hole models. LQC permits 
the exploration of the effects of quantum physics and quantum geometry in 
gravitation p!0| ITT]. 

In order to test the robustness of the LQC it is necessary to apply the 
methodology to some concrete situations and one of the most favorable model 
is represented by the simplest of anisotropic models, namely BI cosmologies. 
The detailed formulation for LQC in the BI models [121 [131 [HI [151 [16] reveals 
the fact that gravity can behave repulsively at Planckian energy densities 
leading to the replacement of the big bang singularity with a big bounce. 

The plan of the paper is as follows: In Section 2 we outline the classical 
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equations for a BI string cosmological model in the presence of a magnetic 
field. In Section 3 we describe the effective loop dynamics for the present BI 
cosmological model. In the next section we present some numerical simula- 
tions and compare the classical and LQC approaches. Some conclusions and 
open problems are discussed in the last Section. 



2 Classical equations 

The line element of a BI Universe is 

ds^ = —dt^ + a\dx^ + a\dy^ + a^dz^ , (1) 

with three scale factors (i = 1, 2, 3) which are functions of time t only and 
consequently three expansion rates. In principle all these scale factors could 
be different and it is useful to express the mean expansion rate in terms of 
the average Hubble rate: 

3^ 3Vai as 03/ 3V' ^' 

where we have defined a new function 

V = = ai«2a3 , (3) 

which is in fact the volume scale of the BI space-time. Hi are the so-called 
directional Hubble parameters: 

= ^. (4) 

In ([2]), dl]) and further over-dot means differentiation with respect to t. 
The Einstein's gravitational field equation has the form 

0,2 03 (32 da 1 , , 
\ \ = -kT-^ , (5a) 

02 03 02 03 

h + '^+h^ = _^T2^ (5b) 

03 Oi 03 tti 

^ + ^ + ^^ = -KTi, (5c) 

Ctl O2 ^1 ^2 



Oi 02 ^ 02 % , % fll _mO 



+ —— + —— = -kT^, (5d) 

Oi 02 O2 O3 O3 Oi 
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where n is the gravitational constant. 

The energy momentum tensor for a system of cosmic strings and magnetic 
field in a comoving coordinate is given by 

= PstringU^u'^ — Xx^X^ + E^^ , (6) 

where Pstring is the rest energy density of strings with massive particles at- 
tached to them and can be expressed as Pstring = Pp + \ where pp is the rest 
energy density of the particles attached to the strings and A is the tension 
density of the system of strings [6l [T71 [18] which may be positive or negative. 
Here Ui is the four velocity and Xj is the direction of the string, obeying the 
relations 

UiU^ = —XiX^ = —1, UiX^ = . (7) 

In (E]) Efj^iy is the electromagnetic field given by Lichnerowich [T^. In 
our case the electromagnetic field tensor F^^, has only one non-vanishing 
component, namely 

F23 = h, (8) 

where h is assumed to be constant. For the electromagnetic field one gets 
the following non-trivial components 



1 B'^ 

p]0 - - -E"^ - -E^ - — = -— 



where /2 is a constant characteristic of the medium and called the magnetic 
permeability. Typically fi differs from unity only by a few parts in 10^ (/2 > 1 
for paramagnetic substances and ft < 1 for diamagnetic) . 

Choosing the string along x^ direction and using comoving coordinates 
we have the following components of energy momentum tensor |20j : 



Tq + pstring — + X — -T^ — -T^ — ■ (10) 

Taking into account the conservation of the energy-momentum tensor, 
i.e., T^.^ = 0, after a little manipulation of (fTOjl one obtains [HI 122] : 

Pstring ~\~ TJ Pstring X . (H) 

V (X\ 

Here we take into account that the conservation law for magnetic field fulfills 
identically. 
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It is customary to assume a relation between Pstring and A in accordance 
with the state equations for strings. The simplest one is a proportionality 
relation [6]: 

Pstring = oA . (12) 

The most usual choices of the constant a are 

1 geometric string 

Q!='^ 1 + u CO > 0, p string or Takabayasi string (13) 
— 1 Reddy string. 



From eq. (]TT|) with (]T2|) we get 

Pstring ~ RO'i'^ 0-2 '^3 ' (-^'^) 

with R a constant of integration. 

Let us also write the other features of BI metric such as expansion and 
shear. The expansion for the BI metric takes the form 

^=^ + ^ + ^ = ^, (15) 

fll ^2 Ct3 V 

while the nonzero components for the shear tensor read 

a, ^ a] = ^ - U. (16) 

In flTBl) and henceforth there is no summation over repeated index "i". The 
shear energy density in given by 

= \a,,a^^ = i {{H, - H,f + {H, - H,f + {H, - H,f) . (17) 



3 Effective loop quantum dynamics 

In the loop quantum cosmology approach we shall use a Hamiltonian frame- 
work where the degrees of freedom of the Bianchi type-I model are encoded 
in the triad components pi and momentum components Cj as follows: 

Pi = 0203, P2 = oiaa, P3 = aia2, Ci = -ydi . (18) 
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Here 7 is the Barbero-Immirzi parameter and represents a quantum ambigu- 
ity of loop quantum gravity which is a non-negative real valued parameter. 
In terms of these variables, the total Hamiltonian of the model is 



'J~C "^grav ~l~ •^matter 

-1 



{C2C3P2P3 + C1C3P1P3 + C1C2P1P2) + y/PlP2P3pM , (19) 



where pm is the matter energy density. In our model pm comprises the 
contribution of cosmological string density Pstring given by ( |T4l) : 



Pstring = RPl (^2^3)^, (20) 

and the energy density of the magnetic field ([9]) [IE] 



P"^^' ~ 2 (a^as)^ 2 pi ' ^''^ 

Einstein's equations are derived from Hamilton's equations: 

d'Ji. d'Ji, 
Pi = {pi, :K} = -/t7^i Ci = {ci, :K} = K'j— . (22) 
dci dpi 

On the other hand, the total Hamiltonian 'K is of constrained type whereby 
it vanishes identically for solutions of Einstein's equations 

:k = . (23) 

Using the explicit form of the Hamiltonian we have for pi the following 
equation 

^ = ^^{cm + C3P3). (24) 

dt 7y/PlP2P3 

and similar equations for p2 and P3. For the evolution of q we get: 

dci Ci 



dt ly/PlP2P3 



{C2P2 + C3P3) 



^lPl^/PlP2P3 



K7 
Pi 



2a V Pi / ^ \pI ) 



(25a) 
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dC2 

Itt 



+ 



C2 



1 



(CiPi + C3P3) 

(C2C3P2P3 + C1C3P1P3 + C1C2P1P2) 



+ — 



K7 
P2 



(25b) 



(it 



C3 



+ 



1 



(CiPi + C2P2) 

(C2C3P2P3 + C1C3P1P3 + C1C2P1P2) 



+ — 



K7 



1 2 ^P2^Y 



(25c) 



Let us observe that from equations for pi and Cj we have the following 
relation: 



^(PiCi) = K.ly/PlP2P3 (^PM +Pi J 



dpM \ 



The directional Hubble rates now reads 



Hi 



PiCi 



ly/PjPk 



i=^j^k = l,2,3. 



(26) 



(27) 



Hamiltonian f|T9l) on account of the vanishing condition (!23|) and fl27|) leads 
to 

i/ii/2 + H^H■, + i/2i/3 = '^Pm, (28) 

which is in fact the zero- zero component (I5d[) of the Einstein system of equa- 
tions dS]). 

Taking into account the symmetry of the density pm with respect to 
variables p2 and ^3 we have 



^(P2C2 -P3C3) = 0. 
at 

This means that for the directional Hubble parameters H2 , H3 we have 

^23 "23 



Ho — H9. 



01^203 A/P1P2P3 ' 



(29) 



(30) 



with a constant. 

The quantum effects in loop quantum cosmology arise in the effective 
Hamiltonian constructed from the classical one by replacing the classical q 
terms with sine functions: 

, sin(^^Q) 

Ci . (31) 

where fli are real valued functions of the triad coefficients pi. 
The effective Hamiltonian is given by: 



^eff = ^ ^ - 1 77~> ^2P3 + cyclic terms |> + ^/piP2P^pM , 



^1 / sin(^'2C2)sin (^^C3) 

(32) 

and the Hamilton's equations for pi and Cj will be modified accordingly. 

It is quite evident that in the limit /ij — )■ 0, the classical Hamiltonian D-C 
(ITU]) is recovered. The expression of the parameters p.[ as functions of the 
triad components Pi represent an ambiguity of the quantization. Two most 
preferable constructions are discussed in [131 El 115] . 

In what follows we shall adopt the /x'-scheme in which the parameters fi'^ 
are chosen as follows: 



with A a constant related to the minimum area gap in LQG. For the numer- 
ical simulations it is assumed that A = 0(1) in Planck units. 

Let us remark that from the vanishing of the Hamiltonian (123]) we have 
the bound: 

^ 1 \ P2P3 . P1P3 . pmX 

P1P2P3PM < — o i — + — + — ( ■ (34) 

I^T I /^2/"3 fJ'lfJ'3 P1P2 J 

In particular, in the fi' scheme the total density is bounded by the critical 
value: 

pMcrii = 3(k7^A)"^ (35) 

implying that the classical singularity is never approached. Indeed the total 
energy of the matter must be below this value and the classical collapse is 
replaced by a bounce. 
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4 Numerical simulations and a comparison of 
the approaches 

The complexity of the equations does not allow for analytic solutions and 
imposes numerical simulations. The classical equations of motion given by 
(I21|) - fl25]) and the corresponding ones for quantum effects in LQC approach 
with the replacement (!3T]) can be solved once the initial values pi(t = to) and 
Ci{t = to) are given. 

In what follows we report some numerical studies on the behavior of pm, 
V and the anisotropic shear cr^,^. Taking into account that only diagonal 
components of shear tensor are non-zero, in the new variables they now read 

= ^ 1 ,^,^t^l,2,3, (36) 

and the shear energy density is 

= [(piCi - P2C2)^ + (P2C2 - PsCs)^ + (P3C3 - PlCi)^] . (37) 

6YP1P2P3 

In doing so we considered a number of cases that helps us to clarify the 
role of various parameters. To begin with we examined both positive and 
negative a. In particular we considered the case with a = 2 and a = —2 and 
it was found that the value or sign of a leaves the overall picture qualitatively 
unchanged. As a second consideration we set a large value of R and small 
value of /? and vice versa, for example, R = 18.24, /? = 1 and R = 0.90, (3 = 
6.5, respectively . It was established that for both cases the overall picture 
remains qualitatively unaltered. Finally we consider the case setting different 
initial values for Cj = Cq. 

In Figs. [Il[3l[5l[7]we have illustrated the evolution of volume scale V (red 
solid line), energy density pm (blue dash line) and shear energy (black 
dot line), whereas in Figs. [21 HJ [6l [8] we have plotted the evolution of the 
components of shear tenor ai (red solid line), a2 (blue dash line) and 
(black dot line). 

In the classical case the initial condition cq < leads to a collapsing 
Universe [cf. Figs. [H [2], while cq > gives rise to an expanding one [cf. 
Figs. El H. 

In the LQC approach, even for Cq < the energy density remains bounded 
below the critical energy (155]) as expected from analytical considerations [cf. 
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Fig. After the bounce, the regime is an expanding one and the Universe 
isotropizes [cf . Fig. Ej . On the other hand, for cq > we have only expanding 
phase of the Universe, the shear energy density TP' remains finite [cf. Fig. [7] 
and again ai tend to zero [cf. Fig. [8]. 

5 Conclusions 

In this report within the scope of Bianchi type-I cosmological model we 
investigate the role of cosmic string and magnetic field on the evolution of 
the Universe. In doing so we employ both classical and LQC approaches. 
It is found that the qualitative picture of evolution does not depend on the 
cosmic string (a), though the value of a leads to quantitative changes. On 
the other hand magnetic field together with cosmic string (given by the pair 
/3 and R) also leave the qualitative feature unaltered. Only initial value of 
Cj's, i.e., initial rate of change of the metric functions a^'s play essential role. 

In the classical approach the initial condition cq < corresponds to a 
classically collapsing Universe, while Cq > is associated with an expansion. 
In the LQC approach for cq < the singularity is avoided via a bounce. 
After the bounce the Universe enters an expansion phase with an asymptotic 
isotropization. For positive cq we get always expansion and isotropization. 

At the classical level other Bianchi models have richer phenomenology 
than the cosmological BI model. From this point of view the extension of 
the string cosmological model in the presence of electromagnetic fields to 
other types of anisotropics and comparisons between the classical approach 
and LQC one deserve further studies. 
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Figure 1: Classical: Evolution 
of pm and Y? for a = —2, 
R = 18.24, (3 = 1, Co = -1. 
Here and further the red solid 
line corresponds to V, blue 
dash line to Pm and black dot 
line to S^. 
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Figure 2: Classical: Evolution 
of fJi's for a= -2, R= 18.24, 
/3 = 1, Co = —1. Here and 
further the red solid line cor- 
responds to C7i, blue dash line 
to C72 and black dot line to C73. 
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Figure 3: Classical: Evolution 
of V, Pm and for a = —2, 
R = 18.24, /3 = 1, Co = 1. 
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Figure 4: Classical: Evolution 
of (Ti's for a = -2, R= 18.24, 
(3 = 1, Co = 1. 
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Figure 5: LQC: Evolution of 
V, pm and for a = —2, R = 
6.20, /3 = 1, Co = -1. 
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Figure 6: LQC: Evolution of 
cr.'s for a = -2,R = 6.20, /3 = 
1, Co = -1. 
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Figure 7: LQC: Evolution of 
V, Pm and for a = —2, i? = 
i? = 6.20, /? = 1, co = 1. 




Figure 8: LQC: Evolution of 
a,'s for a = -2, i? = 6.20, /3 = 
1, Co = 1. 
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